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ABSTRACT 



Using the recent relatively precise experimental results on the pion structure function, obtained 
from Drell-Yan processes, we quantitatively test an old model where the structure function of 
any hadron is determined by that of its constituent quarks. In this model the pion structure 
function can be predicted from the known nucleon structure function. We find that the data 
support the model, at least as a good first approximation. 
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1 Introduction 



In 1974 a model was proposed (T[| for the deep inelastic scattering structure functions of a 
hadron in terms of constituent quarks with structure. For example, the proton is described in 
terms of three UUD constituents with an SU(6) inspired wave function. At large Q 2 the virtual 
photon probes deep into one constituent and sees its parton structure. The proton structure 
function is obtained as a convolution of the Q 2 independent constituent wave function with 
the Q 2 dependent constituent structure function. Similar models of the nucleon in terms of 
constituents with structure have been considered over the years also to describe the static 
properties of nucleons @, ||, In our case, the nucleon structure function is treated in full 
analogy with the case of Helium 3, with constituent quarks replacing nucleons. Of course one 
may object that nucleons, i.e. the constituents of Helium-3, are colourless and therefore can 
exist as unconfmed units. On the contrary, the constituent quarks are confined, so that they 
cannot be really independent of each other and a colour field string must connect them to each 
other. However, it is conceivable that the content of the string in terms of sea partons and 
gluons could be small in comparison with the structure of the constituent. Alternatively, a 
string segment could be associated with the constituent in a universal way, independent of the 
constituent flavour and of the hadron, so that, in a sense, it becomes a part of the constituent 
itself. At the other extreme, the string could be responsible for the whole structure of the 
hadron. In this extreme case we would have a model of three structureless valence quarks and 
a sea of quark and gluons from the string || in principle different for different hadrons. The 
real hadron will probably be somewhat in between. 

In this note we discuss the quantitative information that can be obtained on this issue from 
the available data on the pion structure function which have been collected from measurements 
of the Drell-Yan lepton pair production cross-section. In the proposed model where all the 
structure is in the constituents, one can start from the known parton densities in the nucleon, 
deconvolute the wave function and obtain the parton densities in the constituents. From these 
one can then predict the pion structure function, given a reasonable wave function for the pion. 
We will compare the predictions of this model with the data on the pion structure function 
obtained from the Drell-Yan process ||. This kind of comparison has already been done long 
ago [§|, |7| but the data on the pion structure function || are by now sufficiently precise to make 
the present re-analysis worthwhile. In fact the validity of the constituent-with-structure ansatz 
can now be significantly tested. We shall see that the model is in reasonable albeit not perfect 
agreement with the data. 

It is true that the choice of the wave functions introduces some ambiguity in the prediction 
of the pion structure function from that of the proton. But there are important sum rules in this 
model that are valid independently of the wave function form. In fact the amount of momentum 
carried by gluons, by sea and by valence should be separately the same in the nucleon and in 
the pion at the same Q 2 (as a consequence of the fact that the constituents carry the totality 
of the hadron momentum). No such equality is predicted by the model where all the sea and 
gluons, or a substantial part of them, arises from the string, the structure of the string in the 
proton and in the pion being in principle different. The separate determination of sea and 
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gluons in the pion is difficult, because the available Drell-Yan data do not give any information 
on the pion structure functions at x < 0.2. Most of the information on the gluon distribution 
in the pion arises from the limited data on large px photons produced in n + p reactions |§ . But 
the total momentum carried by sea and gluons is well determined being the complement to 1 
of that of valence and one finds 0.61± 0.02 for the proton |§ and 0.54±0.04 for the pion ||, at 
Ql = 4 GeV 2 . The results of this analysis appear to support to a fair degree of accuracy the 
const it uent s- wit h- struct ure mo del . 



We recall that another application of the formalism of Ref. |]J is for nuclei. In Ref. UTOfl , it 
is shown that a substantial part of the EMC effect (the A dependence of the nucleon structure 
functions) can be attributed to the distortion of the constituent wave function inside a nucleon 
due to the external nuclear field. The model can also be applied to polarized deep inelastic 



scattering [11]]. In this model the constituents carry the totality of the proton spin and the 
observed spin crisis is described by a corresponding depletion of the fraction of the constituent 
spin which is carried by parton quarks. In this picture it is particularly clear that the experi- 
mental results are not at variance with the constituent model. Rather they have implications 
on the constituent structure. 



2 The model 

For definiteness consider a proton p or a positively charged pion tt + . In the model where the 
structure of the hadron is due to the structure of the constituents, the parton density r h (x, Q 2 ) 
for a given parton type r in the hadron h is given by H: 



rh(xiQ ) 



1 dz 

z 



U h (z) ru (- z ,Q 2 ^)+D h (z)r D g,Q 2 ) 



(1) 



where Uh is the density of up-constituents in the hadron h = p,n + , D p is the density of down- 
constituents in the proton p, while D n + is the density of D (antidown) constituents in the pion 
7r+ ) r c/,D are the parton densities in the U or D constituents (for the pion td is actually rg> 
which is the same as fjj). As the Q 2 evolution matrix does not depend on the target, i.e. it 
is the same for the partons in a proton or in a constituent, it follows that if the convolution is 
valid at one Q 2 it will remain valid at all Q 2 . Note that the moments 

r ( :\Q 2 )= I' dxx n -\ h {x,Q 2 ) (2) 
Jo 

are simply given by a sum of products of moments: 

r ( :\Q 2 ) = uPr^{Q 2 ) + D h n) r^{Q 2 ) . (3) 
For short hand we indicate the above convolution by 

r h = [U h ®r v + D h ® r D \ . (4) 
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For example, the gluon density in the proton is given by 

9 P = [U p (g)gu + D p (g) g D ] = (U p + D p ) ® g v (5) 

where the last step is due to the equality of the gluon density in U and D constituents. Actually 
it is important to note that, by using obvious isospin relations (like ur> = du, etc.), for all 
partons kinds r p one can refer to the densities of partons in the U constituent. This is also true 
in the pion case. Now recall that the first moment of U and D are Ujp = 2, = 1 while, 
for the second moments, + D p 2 ^ = 1, because constituents carry the totality of charge and 

momentum of the proton. Hence gj® = g^ . Clearly, for similar reasons, also gf® = gjp . Thus, 
independent of the wave functions, the total momentum of gluons in p and in 7r + are predicted 
to be the same at the same Q 2 . By an identical argument, the same prediction holds for the 
total sea second moment and consequently for the total momentum carried by valence. 



3 Parameters of the model 

In order to predict the parton densities in the pion from those in the proton, we take the proton 
parton densities given by the most recent fits of all available data obtained by Martin et al. in 
Ref. 0. Precisely we use the set of parton densities labeled by MRS(G), with the two-loop 
Q 2 evolution evaluated for A = 255 MeV, where A refers to Nf = 4 in the MS definition, 
corresponding to a s {mz) = 0.114. Other available sets of structure functions will be used to 
check the stability of the results (also with different values of A). As for the distributions of 
the U and D constituents in the proton we take those given in Ref. These constituents 
densities, based on SU(6)w® 0(3) at p z — > oo, are complicated and will not be reproduced here 
(it suffices to say that the parameters introduced in Ref. |jj are fixed to the values (3 = 0.44, 
a 2 = 0.8). Simpler choices would also work and we have checked that no important changes 
in the final results are obtained with different starting wave functions. Then for the parton 
densities in the U constituent we adopt the following parametrisation at Qq = 4 GeV 2 



Qu(x,Ql 



9u(x,Ql 



a 



C a 



x (l+d s ) 



+ 5. 



uq 



X 



(6) 



;i - xy D *-v 



Where B(x,y) is the Euler beta function. In the quark formula there is a valence term, only 
present for the u parton quark, and a universal sea term. The differences in the sea composition 
(strange vs. non-strange, u vs. d etc.) are irrelevant here and have been neglected. At small 
x, a stronger behaviour than 1/x for sea and gluon densities, parametrised by the positive 
coefficients d s and d g , has been allowed according to the results obtained by ZEUS and HI at 
HERA. Of course the momentum sum rule imposes a relation among the parameters. 
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Figure 1: Proton parton densities at Ql = 4 GeV 2 . Solid lines refer to the fit of MRS(G) ||; 
dashed lines refer to the present model. The differences are well within the uncertainties of the 
MRS(G) fit. 

The parameters appearing in the above formulae were fitted to reproduce the input parton 
densities in the proton, given the chosen wave function. The values of the parameters obtained 
from the fit are: A = 0.776, C s = 0.5, D s = 3.3, d s = 0.085, D g = 1.3, 5 g = 0.45. The 
comparison between the input parton distributions of Ref. at and the results of the 
fitted distributions for the model are shown in Fig. 1. As seen, given the accuracy to which the 
parton densities are known, a very good fit is obtained. Thus there is no doubt that the proton 
data are nicely consistent with the model. 

For the total momentum fraction carried by valence in the proton at Qq one has: 

V; (2) (Qo) = 0.39 ±0.02 (7) 

where V = [u — u + d — d\. The error has been estimated by reevaluating the moment starting 
from the available recent compilations of the proton parton densities, as shown in Table |T]. We 
also varied the value of A in a range corresponding to 0.110 < a s (mz) < 0.125 using the recent 
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Table 1: Second moment of valence, and gluon, 
some of the most recent parton density parametrizations. 



in the proton at = 4 GeV 2 for 



results of Ref. 
from varying a 



12| . The total error is a combination of the uncertainty at fixed a s with that 
. The difference in using either the input nucleon densities or the fitted 



densities in the U constituent is completely negligible given the quoted error. 



4 Results 



Having derived the parton densities in the U constituent we now proceed to predict the pion 
structure functions. The only ingredient which is still needed is the distribution of U and D 
constituents in the pion. Following Ref. we take 

U w+ = 

(8) 




V lx +{x) = \ — -exp[— 2f3ln 



1 — x 



where the parameter (3, is fixed to the value (3=0.1 in such a way as to approximately have 
x ■ v n (x, Ql) ~ (1 — x) as x — > 1, according to the Drell-Yan-West relation [13|], and v n is the 
valence quark distribution in the pion. 

The predictions for the parton densities in the pion at Ql are simply obtained by convoluting 
according to eqs. [j], the constituent distributions in eq. 8 with the parton densities in the 
constituent specified in eqs. 6 as a result of the proton fit. In particular, as already mentioned, 
the predicted second moment of valence coincides with the result for the proton and is given 
in eq. 0. 

The above predictions should now be confronted with the experimental data from Drell-Yan 
processes 0. The way to extract the pion structure functions from the Drell-Yan data has been 
recently discussed in Ref. ||. As it is well known the Drell-Yan cross-section is obtained by 
a convolution of the parton densities in the proton times those in the pion times the partonic 
cross-section. The latter includes the QCD correction which leads to a quite substantial K 
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factor [14] at the relevant dimuon mass scale (typically between the J ftp and the T). Thus in 
principle in order to extract the pion densities one has to compute the K factor. The authors 
of Ref. H chose to write the K function in the form K(xp,Q 2 ) = K^'K', where is the 
simple K factor computed at one loop accuracy, while K' includes the effect of higher orders 
(including the correction for a possible bad choice of a s in the leading term). In Ref. || K' was 
fitted from the data. This procedure is only justified if K' is really a constant in xf- Indeed the 
fit is sensitive to a constant rescaling of V. In fact the valence is the dominant contribution at 
the rather large values of x where the data for the pion structure function exist and the overall 
scale of valence is normalized by its first moment = 2, where V w = u + d. A consistency 
check is that K', arising from higher orders, should come out reasonably close to 1. In the fits 
of Refs. H K' ends up in a range between 1.1 -f- 1.3. 

If indeed K' is with good approximation a constant, we can take the results of the fits in 
Ref. H as a compact description of the actual data. Recently an almost complete calculation of 
the rapidity dependence of the K factor at two loop accuracy has been performed in Ref. |15| . 
This calculation is not complete because the effect of soft gluon contributions (within a specified 
definition) is not included. We have repeated the procedure of Ref. §| with the available two- 
loop QCD corrections to the Drell-Yan xp differential cross-section. We found that with a 
very good accuracy K' is indeed a constant over the rapidity range of the experiment. Thus 
we can validate the procedure of Ref. ||. Of course further uncertainties on the result of the 
fit beyond the statistical accuracy arise from the rudimentary parametrization adopted for the 
pion densities and from the assumptions made for the sea and gluon densities which the data 
do not much constrain. 

In Fig. 2 we present a comparison between the model fit of the data and the best fit 
obtained in Ref. ||. The theoretical predictions are presented with and without inclusion of 
the K' factor. We see that the model fits the data quite well but at the price of a somewhat 
larger K' . While in the model independent fit K' ranges between 1.1 and 1.3, in the model one 
needs a larger K', K' = 1.3 4 1.6. 

Clearly the model independent fit has a slightly better \ 2 than the model. Also, in the 
model, the resulting values of the K' factor are a bit too large to be really satisfactory. This 
difference of K' factors is a consequence of the different values for the second moment of valence 
found in the model independent fit and for that implied by the model: 

Vj 2) (A GeV 2 ) = 0.46 ±0.04 (fit) (9) 
yj 2) (4 GeV 2 ) = 0.39 ±0.02 (model) (10) 

where, of course, the model value is the same as in the proton. The error attributed to the fit 
is our estimate which takes into account the error within the procedure of Ref. ||, as given by 
the authors, plus the ambiguities related to the assumptions made in the procedure (mainly 
from the parametrisation choice for the pion, the xf independence of K', the value of a s etc.). 
A plot of the resulting structure function of the pion, in the present model, is shown in Fig. 3, 
where it is compared with the fit of M. 
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Figure 2: A comparison of a sample of experimental Drell-Yan data (tt~W reactions) with the 
fit in Ref. || (solid lines) and with the present model (dashed lines). In each case we present 
result with (upper curves) and without (lower curves) the K' factor. Thus both sets of curves 
can fit the data rather well but the model needs a larger K' factor. 



5 Conclusions 



In conclusion, the constituent-with-structure model is shown to provide a reasonably accurate 
description of the pion structure functions as determined by experiment. This is particularly 
remarkable in that the pion is a very peculiar hadron with mass that vanishes in the chiral 
limit. Thus there is a strong indication that the model can actually provide a reasonable first 
approximation of the structure functions of any other hadron for which no data exist. 

Finally, we recall that the second and third moments of the valence parton densities in 
the pion have been estimated in lattice QCD in the quenched approximation [T^|. There 
the result for the second moment was V; (2) (49 GeV 2 ) = 0.46 ± 0.07 which corresponds to 
VP{A GeV 2 ) = 0.55 ± 0.08. This is a rather large value in comparison with the prediction 
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Figure 3: The resulting valence distribution x ■ v n (x) in the model (solid line) compared with 
the fit of Ref. || given as a band that includes the uncertainties on their procedure. 

of the model. However, it is known that quenched lattice calculations fail to reproduce the 



momentum fractions carried by up and down quarks in the proton (T^|. For the proton the 
lattice results are larger then the fitted values by at least a factor of two. Thus the conclusion 
is that the quenched approximation appears to be rather poor for the calculation of hadronic 
structure functions. 

We warmly thank W.J. Stirling for providing us the fortran code of the MRS(A') and 
MRS(G) parton density parametrizations, P.J. Rijken and W.L. van Neerven for the two loop 
Drell-Yan cross section program and G. Martinelli for the two loop Q 2 evolution program. 
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